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INTRODUCTION 

The  problem  we  are  addressing  In  this  report  Is  how  we  can  design 
signals  which  have  small  volume  contributions  near  the  origin  and  in  a 
strip  parallel  to  the  time  axis  of  the  ambiguity  function.  Because  of 
the  difficulties  of  direct  synthesis  of  signals  from  the  signal  ambiguity 
function,  the  usual  approach  is  to  select  an  interesting  set  of  realizable 
waveforms  and  investigate  the  properties  of  the  selected  set  of  waveforms. 
Here  we  investigate  signals  which  are  composed  of  several  components  of 
the  form 

N 

u(t)  -  u.j,(t)  • 

K*  1 

The  rationale  is  that  if  we  can  choose  a  set  u  (t)  whose  crossambiguity 
functions  appear  in  regions  of  no  interest,  then  this  will  have  the  effect 
of  maximizing  the  signal-to- interference  ratio  because  we  minimize  the 
overlap  of  total  ambiguity  volume  and  clutter  scattering  function. ^ 

The  scattering  function  of  interest  is  a  strip  parallel  to  the  time  axis. 

All  practical  transmitters  operate  under  peak  power  constraints.  In 
order  that  we  can  transmit  the  maximum  possible  amount  of  energy  in  avail¬ 
able  transmission  time,  we  must  design  signals  which  have  the  maximum 
average  power  for  fixed  peak  power.  In  the  next  section,  we  define  the 
concept  of  signal  efficiency  which  will  be  used  to  rule  out  certain 
classes  of  multicomponent  signals  from  consideration.  The  problem  here  is 
that  although  we  gain  clear  regions  in  the  ambiguity  plane,  we  must  reduce 
average  power  (i.e.,  signal  energy)  and,  hence,  little  improvement  is 
realized  in  the  signal  to  interference  ratio  when  one  is  partially  noise 
limited. 

Finally,  we  will  discuss  a  set  of  signals  which  have  100%  efficiency 
and  which  can  be  used  to  move  the  ambiguity  volume  out  of  the  strip  of 
interest.  These  are  the  quadratic  congruential  codes.  The  quadratic 
congruential  codes  are  based  on  ideas  that  have  been  used  in  development 
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II. 


(2) 

of  linear  congruenCial  codes.  We  will  use  the  number  theoretic 

(3) 

properties  of  congruences  to  establish  bounds  on  the  ambiguity  function 

in  all  regions  of  the  time-frequency  plane.  The  results  show  that  as  we 

increase  the  available  time-bandwidth  product,  we  can  move  more  ambiguity 

volume  outside  the  interference  strip  and,  hence,  increase  the  signal  to 

interference  ratio.  Bounds  on  the  volume  and  height  distributions  of  the 

(4) 

ambiguity  function  have  been  established  by  Price  and  Hofstetter. 

What  we  have  developed  in  this  paper  are  easily  generatable,  unity 
efficiency  codes  that  move  the  signal  ambiguity- function  volume  away 
from  the  strip  parallel  to  the  time  axis  of  the  ambiguity  function  (i.e., 
we  eliminate  the  effects  of  sources  of  interference  in  frequency  band 
around  desired  frequency).  This  concept  has  application  to  clutter 
rejection  and  spread  spectrum  communication. 

EFFICIENCY 

The  design  of  the  energy  efficient  codes  has  been  addressed  by 
(7)  (8  9) 

Schroeder  '  and  Ackroyd.  ’  ’  We  define  the  signal  efficiency  as 

Ef  ’  Pavt/Ppeak  *  100%  (2) 


for  a  signal 


u(t)  i  0;  0  <  t  <  T. 


Thus,  we  consider  time  limited  signals.  We  have  that 


P 

ave 


1 

T 


T 

r 


o 


|u(t)  |  2dt 


and 


peak 


MAX  , 
t  |u(t)| 


M 


(3) 


(4) 
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We  assume  the  signal,  u(t),  is  energy  normalized. 


T 

r 


E  = 
u 


i  u (t) I zdt  =  1. 


0 


Thus, 


Hence, 


P 

ave 


1/T 


and 


Ef  =  (1/TM2)  1007,  .  (5) 

We  now  consider  two  classes  of  signals. 

(a)  Single  FM  Signal 

These  are  signals  that  have  only  one  frequency  component  at  one  time. 

They  may  have  many  frequency  components  but  these  components  must  appear 
in  different  time  slots.  We  may  characterize  these  signals  with  a 
rectangular  envelope  as 

u(t)  =  ej9(t),  0<t<T.  (6) 

/T 

Thus,  M  =  — -  and  the  efficiency  is 
✓  T 


"ff 


x  100  =  100%  . 


rr 
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(b)  Multiple  FM  Signals 

At  first  consideration,  these  signals  appear  to  have  many  desirable 
properties.  They  do  not  have  the  pedestal  behavior  that  many  of  the 
single  FM  signals  possess,  such  as  SQFM  and  VCHIRP.  They  have  reduced 
volume  in  the  strip  since  each  FM  component  can  be  designed  so  that 
their  crossambiguity  functions  appear  outside  the  strip  of  interest. 

However,  they  have  a  diminished  efficiency.  Consider 

N 

Zj0K(t) 

aKe  *  0  <  t  <  T.  (7) 

K  =  1 

If  we  assume  for  simplicity  that  the  signal  is  energy  normalized,  each 
component  of  the  signal  has  the  same  energy,  and  the  frequency  components 
are  far  removed  from  one  another  then 

E  -  NTa2  *  1 
u 


and 


(8) 


Hence,  the  maximum  value  for  the  signal  is  at  most 

M  -  Na  -  /N/T  .  (9) 

This  occurs  when  all  cosines  go  through  their  peak  value  at  the  same 
time.  Thus,  we  have 


E 


ff 


N/T 


x  100 


100.. 

N  /s 


(10) 


so  that  we  reduce  efficiency  inversely  and  the  number  of  components  we 
have  in  the  signal. 
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III. 


Computer  results  have  shown  that  for  N  ■  2,  with  two  CHIRP  signals 
spaced  apart,  we  do  indeed  see  a  reduced  volume  in  the  strip  by  a  factor 
of  about  2.  However,  since  we  must  reduce  signal  energy  by  about  the 
same  factor,  virtually  no  gain  in  SNR  is  achieved  under  noise  limited 
conditions.  Thus,  to  prevent  possible  loss  in  SIR,  the  efficiency  of 
signal  must  be  100%. 

MULTICOMPONENT-DISJOINT  TIME  SIGNALS 

It  is  clear  from  the  previous  discussion  that  any  multiple  component 
signal  we  consider  must  also  have  a  100%  efficiency.  Thus,  we  conclude 
that  the  component  must  be  disjoint  in  both  time  and  frequency  in  order 
to  achieve  this  end. 

In  this  section  we  shall  show  that,  if  the  components  are  widely 
spaced  in  time  and  frequency,  we  can  move  the  volume  away  from  the  origin 
and  we  will  discuss  how  much  of  the  volume  is  involved.  In  the  next 
section,  we  present  a  new  class  of  signal  which  provides  our  greatest 
hope  of  achieving  maximum  SIR. 

We  first  consider  two  components  and  then  generalize  the  results  to 
N  signal.  Exact  placement  of  the  components  in  time-frequency  space  will 
be  discussed  in  the  next  section. 

For  two  components,  we  have 

V  ( t )  =  u^t)  +  u2(t) . 


The  ambiguity  function  is 


A 


vv 


Aul  +  u2’  U1  +  u2 


=  Au1u1  +  Au^u^  +  Au^  +  Au2u2> 


(ID 


(12) 
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Taking  the  magnitude  squared,  we  have 


I  Avv  |  2  =  UAu^  +  Au2u2)  +  (Au^  +  Au^)]2 

/ - A - v  / - B - \ 

2  2 
=  lAu^  +  Au2u2|  +  lAu^  +  Au2ux| 


f 


+  2Re| 


(Au^  +  Au2u2)(Au1u2 
-  C  - 


+  Au2u^)*  > . 
_ / 


(13) 


We  discuss  the  various  terms  in  Equation  (13)  labelled  A,  B  and  C. 

First,  we  observe  that  since  the  signals  are  widely  spaced,  the 
peaks  in  the  cross  term  and  the  auto  term  appear  in  totally  different 
regions  of  the  time-frequency  plane.  Hence,  C  is  very  small  and  can  be 
ignored.  By  the  same  argument,  B  can  be  shown  to  be 


B  = 


+  |Au2ux 


2 


Thus,  the  total  is 


|AW|2  =  lAu^  +  Au2u2|2  +  |Auxu2|2  +  jAu^J2. 


(14) 


Also,  we  know  that 


| AulU2 (t ,^)  I  dTd<f>  =  Eu2Eu2. 


(15) 


Suppose  for  simplicity  that  each  component  has  the  same  energy.  Thus, 


Eu^  =  Eu2  =1/2  and  Ey  ■  1. 
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Since  the  total  volume  is  unity  and  the  volume  under  each  cross  term  in 
(14)  is  1/4,  we  have  moved  half  of  the  volume  away  from  the  main  lobe  and 
placed  it  at  the  peak  of  the  two  cross  terms. 

Consider  now  an  N-component  signal. 


v(t) 


£  -.<» 


(16) 


Each  component  will  have  Eu^  =  1/N.  By  the  same  argument,  we  have  that 


vv 


2  =  |AU]Ul+. 


+AununI 


N 


+  X.  iA,y/- 


J',K=  1 
K#  j 
2 

the  second  term  of  which  has  N  -N  terms. 


Since  the  signals  are  assumed  widely  separated  and  if  the  cross 
terms  do  not  overlap,  we  have  that  the  second  term  has  volume. 


V  =  (NZ-N) (— j-)  =  1  -  1/N 
sec  N2 


(17) 


(18) 


and,  hence,  we  have  moved  all  but  1/N  of  the  volume  away  from  the  origin. 
The  difficulty  here  is  that  if  the  volume  is  still  in  the  strip  or  if 
the  cross  terms  overlap,  we  have  lost  much  of  the  gain  for  our  present 
application.  For  example,  the  cross  terms  may  themselves  pile  up  consider¬ 
able  volume  although  away  from  the  main  lobe,  yet  in  the  strip.  In  the 
next  section,  we  will  discuss  a  placement  in  T-W  space  which  assures  that 
these  problems  do  not  occur.  We  shall  use  a  guard-band  in  frequency  to 
account  for  the  problem  of  zero-time  cross  terms. 
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IV.  QUADRATIC  CONGRUENCE  PLACEMENT 

In  this  section,  we  discuss  a  method  of  placing  the  signal  components 
in  time-frequency  space  which  assures  that  we  will  have  the  maximum  SIR 
with  100%  efficiency.  The  notion  of  efficiency  requires  that  we  consider 
essentially  pulse-train  like  signals,  each  of  which  possesses  a  different 
spectral  component. 

The  fundamental  notion  that  we  require  is  as  follows:  the  degree  to 
which  two  signal  components  correlate  with  one  another,  if  they  occupy  the 
same  time  slot,  is  determined  by  their  proximity  in  frequency.  If  they 
occupy  the  same  spectral  location,  they  will  maximally  correlate.  As 
they  separate,  their  correlation  decreases  inversely  as  the  frequency 
difference  between  them. 

This  concept  was  used  to  develop  a  set  of  codes  which  have  mutually 

small  crosscorrelation  properties  for  all  shifts  in  time  and  frequency 

(2) 

axis.  These  are  the  linear  congruence  codes.  The  difficulty  here  is 
that  we  are  not  interested  in  a  set  of  codes  but  only  one  placement. 
Further,  the  autocorrelation  properties  of  these  codes  were  such  that 
significant  side  lobe  could  appear  within  the  strip  of  interest  (i.e., 
for  zero  Doppler). 

We  shall  first  require  that  a  guard-band  be  placed  between  the  top 
of  a  low  component  and  the  bottom  of  the  next  component  up.  This  assures 
that  for  the  present  application  we  have  minimum  strip  volume.  We  call 
this  guard-band  G  Hz.  Each  component  will  be  separated  by  at  least  G  Hz 
from  all  others.  Now  we  call  the  spectral  centroid  of  each  component 
fR,  so  that  the  set  {f^,}  is  what  we  must  choose  in  order  to  define  the 
signal.  The  individual  components  themselves  may  be  chosen  to  suit  the 
desired  behavior  at  and  around  the  main  lobe;  whereas,  the  set  (fK) 
determines  how  large  and  where  are  the  side  lobes  of  the  signal.  Since 
most  of  the  volume  (see  Equation  18)  occurs  in  the  side  lobes,  this  is 
the  cricial  choice. 

We  shall  first  demonstrate  a  poor  choice  of  centroid  set.  Suppose 
we  take,  as  our  set,  a  linearly  increasing  frequency  f^  =  f^  +  KAf 
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K=0,...,N-1.  This  is  shown  in  Figure  1  for  N=5. 


if 


Figure  1.  Linear  Congruential  Code  for  N=5. 

There  are  two  problems  with  such  a  signal.  First,  if  we  shift  the  signal 
one  slot  horizontally  and  one  slot  vertically,  then  4  dots  overlap  one 
another.  Similarly,  if  we  shift  two  slots  in  both  directions,  then  3  dots 
overlap.  This  accounts  for  the  ridge  behavior  that  CHIRP  signals  possess 
in  their  ambiguity  functions.  Secondly,  if  we  shift  horizontally  one  slot, 
then  all  but  one  of  the  dots  are  the  same  distance  apart.  This  accounts 
for  the  poor  side  lobe  structure  that  CHIRP  signals  have  in  their  auto¬ 
correlation  function. 

What  we  desire  is  a  signal  that  has  no  more  than  one  overlap  for 
every  shift  and  a  nonconstant  difference  for  horizontal  shifts.  To  this 
end,  we  now  consider  a  quadratic  of  the  form. 

yK+l  =  tyK+  (K+1>]Mod(N>» 

0< K < N-l 

In  what  follows,  we  shall  assume  that  N  is  the  prime  number.  We  also 
assume  that  y^  is  zero.  However,  the  actual  frequency  will  be  f^,  the 
lowest  spectral  centroid.  Figures  2  and  3  show  several  examples  of  the 
arrays  generated  for  several  values  of  N.  Two  of  these  are  prime;  whereas, 
the  third  is  not.  Notice  the  symmetry  with  respect  to  the  center  for  all 
of  these  arrays.  This  may  first  appear  to  preclude  the  use  of  these  sets. 
However,  we  get  around  this  by  only  using  the  first  (or  second)  half  of 
the  arrays. 
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We  shall  now  present  the  properties  of  the  quadratic  congruences 
defined  in  Equation  (19),  leaving  the  proofs  for  the  Appendix.  We  will 
then  apply  these  results  to  establish  the  results  for  the  ambiguity  function. 

prop.  1)  y (K)  =  Mod(N)  0<  K<N-1 

prop.  2)  Symmetry :  for  any  N,  odd 

y(N-K-l)  =  [y(K)J  Mod(N)  0<K<N-1 

prop.  3)  Horizontal  Shift: 

y (K  +  m)  =  fy(K)  +  Km  +  y(ra)]  Mod (N)  0<  K,m  <N-1 

prop.  4)  Difference  with  Horizontal  and  Vertical  Shift: 

Z(K,m,£)  £  y(K  +  m)-y(K)  +  i 

-  [y(m) +Km  +  f,]Mod(N)  0  <  K,m,«.  <_  N  -  1 

Property  one  establishes  a  closed  form  for  generating  the  arrays 
without  the  recursive  formula.  It  shows  why  we  call  these  quadratic 
congruences.  The  formula  is  quadratic  in  K.  Property  two  is  the  symmetry 
property  we  observed  previously.  Property  four  is  the  main  result.  We 
observe  that  for  any  (m,£)  with  m  ^  0,  the  difference,  2,  is  a  linear 
congruence.  Thus,  for  N  prime  the  differences  must  go  through  a  complete 
residue  sequence.  Mod  N.  Thus,  we  know  the  minimum  distance  (in  frequency) 
for  the  entire  time- frequency  plane.  For  the  ra  -  0  case,  all  the  shifts 
are  l  units  apart.  This  is  another  reason  for  the  guard-band.  The 
property  for  horizontal  shifts  of  the  CHIRP-like  signals  switches  axes 
and  now  is  valid  for  strict  Doppler  shifts  with  no  time  shift. 

Since  the  differences  go  through  a  complete  residue  class,  there  can 
be  only  one  intersection  of  the  two  frequency  patterns.  If  we  were 
extending  the  patterns  (Mod  N),  there  would  be  exactly  one  crossing. 

Note  that  for  horizontal  shifts  this  always  occurs  on  the  opposite  side 
from  the  center  (see  Figures  2  and  3),  as  long  as  N  is  a  prime  number. 
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N  =  37 


Figure  3. 


‘  * 


Quadratic  Congruentlal  Codes  for  N  ■  19  (Prime) , 
N  ■  21  (Non-Prime)  and  N  ■  37  (Prime). 
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V.  BOUNDS  ON  THE  AMBIGUITY  FUNCTION 

We  have  established  the  number  theoretic  properties  for  the  set  of 
frequency  differences.  In  light  of  these  results  and  the  results  for 
linear  congruences,  we  can  place  upper  bounds  on  the  height  of  the 
ambiguity  function  for  the  quadratic  congruence  signals. 

Let  us  assume  that  the  component  signals  are  themselves  cosines  at 
the  frequencies  {f and  that  the  total  bandwidth  is  F  (-NG).  The  actual 
components  may  be  SQFM,  CHIRPS  or  anything  else  for  that  matter.  The 
choice  of  signals  only  affects  the  main  lobe  behavior  and  the  shape  of 
the  non  main  lobes,  but  not  their  position  and  amplitude. 

Under  this  assumption,  we  observe  that  two  frequencies,  f  and  f  , 

m  n 

which  overlap  for  T^  seconds  have  maximum  crosscorrelation 


C 

mn 


N/7fT1  |m  -  n|F, 

1/N , 


m^n 


m  -  n 


(20) 


We  are  also  assuming  that  each  element  of  the  set  {fR}  is  one  of  the 
frequencies 


{fQ  +  KG  ;  0  £  K  _<  N  -  1). 

We  are  simply  assuming  that  the  components  are  the  same  distance  apart. 

Property  four  shows  that  for  the  region  away  from  the  zero  time  slot, 
the  frequency  difference  goes  through  a  complete  residue  class.  Thus,  we 
may  bound  the  height  of  the  ambiguity  function  with  where 

(2) 

This  result  is  proven  in  the  linear  congruence  coding  paper.  In  the 
zero  time  slot,  where  the  frequency  differences  are  all  the  same,  we 


-HI 


(21) 
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bound  the  ambiguity  function  with  A2  where 


(22) 


Note  that  A2  is  somewhat  higher  than  A^  since  an  N2  appears  in  the 
numerator.  The  1/N  term  in  A^  comes  from  the  fact  that  there  will  be 
exactly  one  overlap  (Mod  N) ;  for  the  nonperiodic  case  there  will  be 
at  most  one  overlap.  In  the  light  of  property  two,  if  we  were  to  only 
use  the  first  half  of  the  frequencies,  we  can  drop  this  term  and  bound 
the  zero  Doppler  strip  which  includes  the  autocorrelation  function  by 


A1(l/2) 


4N 

7TFT 


(23) 


This  result  is  valid  away  from  the  main  lobe  (m  t  0) . 

Now  we  can  see  the  trade-offs  between  removed  volume,  number  of  slots, 

and  time-bandwidth  product.  Thus,  as  we  push  more  volume  outside  this 

strip  by  increasing  N,  then  we  must  increase  FT  in  order  to  lower  the 

amplitude  of  the  ambiguity  function.  Thus,  we  have  established  the 

trade-offs  between  strip  volume  and  peak  amplitudes.  The  peaks  away  from 

the  main  lobe  and  off  the  zero  Doppler  and  zero-time  strips  are  governed 

by  A^  which  has  an  additional  1/N  term  in  the  bound.  Thus,  these  peaks 

will  be  slightly  higher  (depending  upon  N) .  Finally,  the  worst  case  occurs 

in  the  zero  time  slot  (where  the  removed  volume  is  going),  indicating  the 

2 

amplitude  is  increasing  as  N  . 


CONCLUSIONS  AND  FURTHER  WORK 

In  this  report,  we  have  established  the  concept  of  efficiency  and 
shown  that  multicomponent  signals  can  be  used  in  clutter  limited 
conditions  only  if  the  components  occur  in  disjoint  time  slots.  We  have 
defined  a  new  class  of  signals  which  provides  a  mechanism  for  removing 
volume  from  a  strip  about  the  zero  Doppler  axis  and  also  placing  bounds 
on  the  actual  heights  of  the  ambiguity  function.  The  trade-offs  between 
volume  removed  and  ambiguity  amplitude  have  been  established  for  this 
class  as  well. 
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VII. 


What  remains  is  to  choose  the  size  of  the  guard-band  and  the  various 
component  signals.  The  guard-band  is  chosen  by  the  strip  size  we  wish  to 
clear.  This  is  determined  by  clutter  data.  The  component  signals  are 
chosen  to  improve  main  lobe  behavior  of  the  ambiguity  function  of  the 
multicomponent  signal. 

We  are  presently  investigating  the  role  of  the  component  signals  and 
the  shapes  of  the  ambiguity  functions  for  various  component  signal  choices. 
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Quadratic  Congruences 
(1)  We  define  y(K)  by 

y  (K  +  1)  «  [y(K)  +  (K+l)]  Mod  N 

y(0)  =0,  0  <  K  <  N  1,  N  prime 

"use  a  as  "is  congruent  to.  Mod  N" 


(2)  T.  y  (K)  a 


_  [~K(K+1) 

=  L  2 


is  a  solution  to  (1) 


P.  Consider  y(K+l)  s  2>J 


.  ptLUj  +  x  iK+u  Q.E.„. 


(3)  T.  For  N  odd  y(K)  s  y(N-l-K)  (symmetry  about  center  point) 


P.  y(N  -  i  -  K) 

,  |~CN-1-K)(N-K)j 

expanding  .  [SflUO*  .  i2S±ii  +  £] 

-  p^N~  ^  -  nk  +  k^k  ~  in 


However,  A  a  0,  B  a  0.  Thus,  since  C  a  y(K) 
y(N-l-K)  =  y  (K) 
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(4)  Translation  in  K 

T.  y(K  +  m)  s  y(K)+Km  +  y(m) 


P.  y(K  +  m)  =  (K  +  m)  (K  +  m  +  1)  g  j~K(K+l  +  m)  +  m  (K  +  m  +  1)  j 

f~K(K+l)  Km  m(tn  +  1)  Km"] 

L  2  2  +  2  +  2j 

,  pl£tU  +  k»  +  aia±ii]  q.e.d. 


(5)  Corr.  Add  £-y(K)  to  both  sides 

Z(K,m,£)  =  y(K  +  m)  +  £-y(K)  =  y(m)  +  Km  +  £ 

/ 

Thus,  for  fixed  translations  (m,£),  m^O,  Z(K,m,£)  is  a 
linear  congruence  in  K.  Hence,  Z  goes  through  a  complete 
residue  sequence  as  0  <_  K  <_  N-  1.  For  m  =  0,  Z(K,£)  =  £. 
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